Kondo Insulator description of spin state transition in FeSb 2 
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The thermal expansion and heat capacity of FeSb2 at ambient pressure agrees with a picture of 
a temperature induced spin state transition within the Fe t2 9 multiplet. However, high pressure 
powder diffraction data show no sign of a structural phase transition up to 7GPa. A bulk modulus 
B=84(3)GPa has been extracted and the temperature dependence of the Griineisen parameter has 
been determined. We discuss here the relevance of a Kondo insulator description for this material. 

PACS numbers: 75.30.Mb, 71.28.+d, 75.20.Hr 
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FcSb2 is a small gap semiconductor whose magnetic 
properties strongly resemble FcSi. The magnetic suscep- 
tibility of FeSi evolves from a high temperature maximum 
to a flat van-Vleck like behavior as T^O with a Curie 
tail which is sample dependant and attributed to defects 
and impurities. Acppli and Fisk have pointed out that 
the underlying physics of FeSi is similar to the physics 
of Kondo insulators with localized magnetic momcntsi 
This picture was also supported by other measurements 
as welli^^ All Kondo insulator materials are cubic with 
the exception of CeRhSb and CeNiSn and are all 41 VCir 
termetallics, except FeSi, RuAl2 and Fe2VAlii4«i A sim- 
ple model that explains their properties involves a flat f 
band hybridizing with a broad conduction band with two 
electrons per unit cell. At T=0 the electrons populate a 
lower hybridized band, and with the increase of T the 
electrons start populating the higher band, resulting in 
a thermally activated Pauli susceptibility. The Anderson 
Kondo lattice model with two electrons per site with one 
/ and one conduction band thus captures the important 
physics. Cubic structures are more favorable for this sce- 
nario since the requirement of just one band crossing the 
Fermi level is more likely. 

Another way to think about FeSi is based on the pic- 
ture of a nearly magnetic semiconductor in an itinerant 
model&£ This picture was supported by local density 
approximation (LDA) band structure calculations^^ An 
important question to answer is if the itinerant picture is 
correct or if the framework of the Anderson model with 
a reduced on-site Coulomb repulsion U better describes 
the physics of FeSi. Other possible model systems con- 
taining 3d metals where these scenarios can be studied 
are therefore highly sought after. 

The small gap semiconductor FcSb2 crystallizes in the 
marcasite-type FeS2 structure (Fig. l(a))^l The basic 
structural unit is made out of Fe ions surrounded by a 
deformed Sb octahedral The Sb octahedra share edges 
along the c axis (Fig 1(b)). According to "classical" 
ligand field theory, the Fe t2 S orbitals are further split 
into two lower lying A orbitals with d yz and d xz sym- 
metries and higher lying S orbitals with d xy symmetry, 



which are directed toward near-neighbor cations along 
the c-axis of the crystal^ One type of marcasite the 
so-called "anomalous marcasite" phase is formed with 
cations having d™ configurations for n>6 and a large c/a 
ratio of 0.73-0.75 in which the d xy orbitals are doubly oc- 
cupied. The other so-called " regular marcasite" structure 
is formed with n<4 for smaller c/ a ratios of 0.53-0.57 be- 



cause the d xy orbitals arc empty 



11-13 



Consequently, the 



Fc-Sb-Fc bond angle a between the neighboring cations 
in edge-sharing octahedra along the c-axis is a < 90° 
in regular and a > 90° in anomalous marcasites (see 
Fig. 1(b)) iiS Thus, if the d n configuration of transition 
metal atom has n > 4, then it implies electrons in H or- 
bitals. Repulsive forces along the c-axis will consequently 
increase the angle a. For example, FeS2 with a d 6 con- 
figuration has a=97.5° as opposed to FeAs2 with a d 4 
configuration and a=72.5°. Therefore the angle a should 
increase monotonically with the number of S electrons 
from a d 4 configuration (S empty) to a d e configuration 
(S filled) ii 

The ionic configuration of iron in FeSb2 is d 4 and its 
ground state at low temperatures is the low spin state 
Fe 4+ (tf g eg) with the A orbital filled and the S orbital 
cmptytiSi Indeed, at low temperature the magnetic sus- 
ceptibility of FcSb2 results entirely from the core dia- 
magnctism and shows little temperature dependence ~ 
Similarly to FeSi, above 100K the susceptibility is 
paramagnetic^^ One possible explanation for the en- 
hanced magnetic susceptibility is a spin state transition 
of the Fe ions^ It is not uncommon that in transition 
metal alloys the crystal field energy marginally exceeds 
the Hund's exchange energy and therefore energies of the 
order of lmeV can produce excitations across a spin gap. 
What makes anomalous marcasites unusual is the exis- 
tence of a band of itinerant electronic states due to the 
overlap of higher lying S orbitals along the c-axis^S In 
this work, structural, thermodynamic and magnetic stud- 
ies were carried out in an effort to further elucidate the 
structure-property relation in this material since a spin 
state transition within the 3c? multiplet might couple to 
the lattice. 
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FIG. 1: (a) Crystal structure of FeSb2 showing Fe (open, 
brown crosses) surrounded by Sb (full blue) octahedra. The 
two short Fe-Sb bond distances are shown in black and white. 
The d xy orbital is oriented in between short Fe-Sb bond dis- 
tances along the z-axis (the c-axis of the unit cell), (b) The 
Fe-Sb-Fe bond angle a associated with the edge sharing oc- 
tahedra along the c-axis of the unit cell. Note that increased 
occupancy of d xy orbitals is reflected in the increase of a. 
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FIG. 2: Temperature dependance of (a) normalized unit cell 
parameters (b) Bond distance lengths and (c) Fe-Sb-Fe angle 
a between edge sharing octahedra and bond distance ratio 
(long/short). Note the increase of angle a above T crossover 
(see text) and the presence of two regions of the Fe-Sb nor- 
malized short bond distance length contraction around it. 
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FIG. 3: Pressure dependance of (a) Normalized unit cell pa- 
rameters (b) Unit cell volume of FeSb2 



High quality FeSb2 single crystals have been grown 
from excess Sb fiux^ Variable temperature and high 
pressure experiments were performed at the beamline 
X7A of the National Synchrotron Light Source at the 
Brookhaven National Laboratory. A closed-cycle He- 
refrigerator and a diamond-anvil cell were used to con- 
trol the sample environment, and monochromatic syn- 
chrotron X-ray and gas-proportional position-sensitive 
detector were used to measure the powder diffraction 
data^& Magnetization and specific heat were measured in 
Quantum Design MPMS and PPMS instruments. Single 
crystals of FeSb2 were thoroughly ground to a fine pow- 
der. Rietveld refinements were performed using GSAS^i. 
The results of fits were summarized in Tables 1 and 2. 

Fig. 2 shows the temperature evolution of sev- 
eral structural parameters. Linear coefficients of ther- 
mal contraction^ a in the paramagnetic stated on = 
l/l(dl/dT)p (i = a, b, c) for individual unit cell axes 
are: a a = 4.2 • lO -6 /*" -1 , a b = 1.3 • lO^A^and a c = 
3.4 ■ 10 _5 A _1 . Individual Fe-Sb bonds also change as 
the temperature is lowered and temperature effects on 
bond distance lengths normalized to 500K are shown in 
Fig. 2(b). We observe two regions of bond contrac- 
tion with a crossover near 100K. The thermal contrac- 
tion of the long Fe-Sb bond distances (in the octahe- 
dral plane) becomes marginally stronger below 300K, 
whereas the short (apical) Fe-Sb bond distances con- 
tract significantly between 500K and 100K. The con- 
traction of the short Fe-Sb bond distances is finished 
at 100K and below this temperature we observe no fur- 
ther thermal contraction on cooling down to the low- 
est measured temperature of 20K. We can estimate the 
bond distance contraction by calculating the mean coef- 
ficient of contraction between 500K and 20K: eETi,T 2 = 
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dicating relatively large differences in the distortion of 
the octahedra. Both the bond distance ratio and the Fe- 
Sb-Fc angle a associated with the shared octahedral-site 
edges along the c-axis (Fig. 2 (c)) change the slope be- 
low 100K. The change in the bond distance ratio implies 
maximum distortions for the Fe-Sb octahedra around the 
^crossover ~100K. The value of the Fe-Sb-Fe bond angle 
at 300K (76.028±0.02)° agrees well with the previously 
reported value of 76. O W 3 Fig. 2(c) shows that the Fe- 
Sb-Fe bond angle a has little or no temperature depen- 
dance for T \ 100K and increases above T cros;!0t , er ~100K. 

The lattice parameters contract smoothly with pres- 
sure to smaller values above 3GPa than the ones ob- 
tained by thermal contraction down to 20K (Fig. 3(a), 
Fig. 4). Individual lattice parameter compressibili- 
ties /3q = —j(-§p) T (i = a,b,c) are quite different: 
/3§ = 0.00223(4)GPa-\ $ = 0.00339(3)GP a - 1 and 
/3g = 0.00655(5)GPa" 1 . The ratio of the c and a axis 
compressibilities is as much as (4^)^=2.9. When com- 
pared with a ratio of thermal contraction (^)p=8, it 
follows that structural changes are more anisotropic with 
temperature. Both temperature and pressure compress 
the c-axis which is the direction of higher lying 5 orbitals 
associated with d xy symmetry rather strongly, whereas 
the b and a axis are less compressible. The ratios of con- 
traction to compressibility are also highest for the c-axis: 

(jfY = LSS-lO^GPaA"" 1 , (|) = S.SS-lO^GPaA"- 1 

and ^\ = 1.15 • 10~ 2 GPaK~ 1 . We see no evidence 

for the formation of a high pressure polymorph up to 
(7.14±0.1)GPa, as was established in the cases of CrSb2 
and NiSb2^ The uncertainty in Rietvcld refinement and 
the errors of the bond distances at high pressures limit 
our conclusions on the bond distance ratio and the Fc-Sb- 
Fe octahedral angle between neighboring cations in edge- 
sharing octahedra. However, preliminary measurements 
indicate that there is no change in the pressure depen- 
dence of the Fe-Sb bond distances and the angle a in con- 
trast to their thermal contraction. Their ratio remains 
constant under pressure. Fig. 3 (b) shows the normal- 
ized unit cell volume as a function of pressure. We obtain 
the value of bulk modulus B =84(3)GPa and B'=5(l) by 
fitting a Birch-Murnagham equation of stated 

The Fe-Sb long and Fe-Sb short bond lengths in Fe-Sb 
octahedra involve Fe a, a bonding orbitals (e g in origin) 
and Sb sp 3 hybrid orbitals, as pointed out by Goode- 
nough and F&nH^ Bond distance lengths in Fe-Sb oc- 
tahedra influence the large energy gap between t2 9 and 
e g derived states. The gap within the t2 g multiplct is 
caused by mixing of S cation and a bonding anion or- 
bitals, and the strength of this mixing is governed by the 
Fe-Sb-Fe angle a that measures the deviation from an 
ideal tetrahedral angled Starting from the diamagnetic 
state state at low temperatures, the increase of the the 
Fe-Sb-Fe angle a between neighboring cations (Fig. 1(b), 
Fig. 2(c)) above T crossover ^100K indicates a thermally 
induced population of higher lying S orbitals (d^j, in ori- 
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FIG. 4: The effects of pressure and temperature on the unit 
cell volume. The 8GPa of pressure reduce the unit cell volume 
by 1.4%, 4.6 times exceeding the ambient pressure compres- 
sion for 6.5% by temperature decrease from 300K to 20K. 



gin). Since T crossover coincides with the temperature of 
diamagnetic to paramagnetic transition in FcSb2 , we con- 
clude that the spin state transition observed in Rcf. 14 
occurs within the Fe t2 9 multipletii^ 

Goodenough's description of the energy bands of FeSb2 
can be connected with the electronic and magnetic prop- 
erties of Kondo insulators. Transitions within the t2 9 
multiplet involves the population of a d^^-derived band 
of itinerant states. In the case of a Kondo insulator we ex- 
pect to describe the magnetic susceptibility with a model 
of a metallic spin paramagnetism, albeit with a large 
renormalization of the noninteracting bands. Thus, we 
used the model proposed by Jaccarino^ The tempera- 
ture dependence of the magnetic susceptibility should be 
explained with a model of two narrow bands (peaks) at 
the density of states of width W and separated by E g =2A 
(see inset in Fig. 5). This model has been successfully 
applied to show the validity of Kondo insulator picture 
in seminal work on the thermodynamics of FeSi£ 

The Pauli susceptibility of an itinerant electron system 
with N(E) density of states is: 



x(T) = -2/4 / N(E) df ^ T) dE 

c—band 

where (is is the Bohr magneton and f(E,fi,T) = 
(cxp[(P - f-i)/k B T] + I) -1 . The factor 2 in the above 
equation is due to the holes in the valence band that 
contribute to \ in the same way as electrons in the con- 
duction bandS Taking N(E) = Np/W and [i = W + A, 
where p is the number of states/cell, and N number of 
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FIG. 5: Magnetic susceptibility of FeSb2 (data taken from 
Ref. 14). Note much smaller low temperature Curie tail than 
in FeSi. Solid and dashed lines represent fits using the model 
density of states shown in the inset. 



unit cells, we obtain: 
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FIG. 6: (a) Heat capacity Cp of FeSb2. Inset shows contri- 
bution to heat capacity from the model density of states used 
to model magnetic susceptibility (see text), (b) Griineisen 
parameter F of FeSb2 
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By fitting magnetic susceptibility data we obtain allow- 
able values for parameters: OK < W/ks < 400K and 
850K < Eg/ks < 1100A', comparable to parameter val- 
ues obtained for FeSid 

The accumulation of states at the gap edges should 
in principle be observed in other thermodynamic prop- 
erties, such as the thermal expansion a and the heat 
capacity Cp. Unfortunately, there is no MSb2 diamag- 
netic member with M=3d transition metal clement with 
the marcasite structure that we could use in order to 
extract the electronic contribution to a and Cp. A simi- 
lar difference has been obtained for FeSi using CoSi as a 
referenced CoSb2 forms in a distorted monoclinic struc- 
ture P2\/c at ambient temperature and NiSb2, contrary 
to previously reported, is not diamagnetici^* From the 
low temperature specific heat (Fig. 6a) C ~ BT 3 where 

B = ^ni?7r 4 (jf^J , we estimate the Debye tempera- 
ture as 9d = (256 ± 6)K, a smaller value than pre- 
viously reported (8d = 380A') based on vapor trans- 
port grown samples. Using the same range of parame- 
ters obtained from the magnetization fit, we have calcu- 
lated spin state transition contribution to specific heat 
Cp=(dU/dT) with U given by the contributions from 
valence and conduction band: 



IV 




Np EdE 
W exp(^(£; -W-A)) + l 



o 

2W+2A 

f Np EdE 
+ J W exp(f3(E - W - A)) + 1 ' 

W+2A 

The results are shown in Fig. 6a where we see that 
phonon contribution dominates at T~T crosso „ er . From 
the measured specific heat (Fig. 6a) and the unit cell 
volume temperature dependance (Fig. 4) we estimate 
the Griineisen parameter T = 3Ba/Cp (Fig. 6b), where 
a = 3^(fp-)p (Fig- 7) and B is the bulk modulus. The 
overall magnitude of the thermodynamic Griineisen pa- 
rameter r is comparable to the one seen in FeSii^ We see 
that r is weakly temperature dependant with a possible 
minimum developing around T=50K. Since T represents 
the measure of anharmonic effects in the lattice potential 
energy, our T(T) indicates the overall agreement with the 
Griineisen assumption Au/w = — TAV/V. 

The aiT) for FcSb 2 (Fig. 7) is different from the case 
of cubic or anisotropic metals, where a(T) follows the 
C p decrease below room temperature in accord with the 
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FIG. 7: Thermal expansion a — y(^)p of the unit cell 
volume of FeSb2. Weak temperature dependance of the unit 
cell volume introduces some noise in thermal expansion data 
below 100K. Inset shows thermal expansion coefficient of the 
individual axes a = j(-§p)p (l=a,b,c) 

Debye function. Departures are expected for large vari- 
ations of The electronic contribution to the specific 
heat of a narrow gap semiconductor is expected to show 
a Schottky peak, as expected for a two level electronic 
system separated by gap E g =2A. If the degeneracies of 
levels separated by E ff are equal, we expect the peak at 
a temperature equivalent to A = 2AT pea k~^ While 
we cannot extract the electronic contribution from the 
phonons in Cp (Fig. 6a shows that it is ~10% of the 
measured value of Cp), since the Griineisen parameter 
r = 2>Ba/Cp remains roughly constant, the Schottky 
term should be reflected in the thermal conductivity as 
well. Notwithstanding that we were not able to extract 
the electronic contribution to a, the data at Fig. 7 could 
be consistent with a broad peak forming at T pea k ~(125- 
200)K, thus giving A/k B ~ (300 - 480)Jf. This might 
indicate that physically meaningful values of the gap are 
850K < Eg/ks < 960-ftT, closer to the lower limit of 
those obtained from the fit of the magnetic susceptibil- 
ity. Anisotropic thermal expansion coefficients of the the 



lattice constants ai = j(-§p)p (I = a,b,c) (Fig. 7 inset) 
show that the peak in the volume thermal expansion is 
dominated by the c-axis contribution a c , whereas we see 
little or no evidence for the peak formation for a a or at,. 
Since the d xy orbitals that could be involved in a spin 
state transition within the t2 g multiplet overlap along 
the c-axis of the crystal^, it is likely that the electronic 
states involved in the Schottky peak formation are in fact 
of d xy origin. This agrees with Goodenough's picture of 
a band of itinerant states within the t2 g multiplet and 
observations of enhanced conductivity in paramagnetic 
region. 12 1 14 1 15 

In conclusion, we have shown that the experimental 
data for FeSbo can be viewed BjS ci thermally induced 
spin state transition within the Fe ti g multiplet. A sim- 
ple model with two peaks at the density of states at the 
gap edges can explain the thermally induced paramag- 
netic moment, as expected for Kondo Insulators. The 
magnetic properties correlate with the electronic con- 
ductivity in FeSb2 since the resistivity in the low spin 
diamagnetic state at 2K is more than four orders of mag- 
nitude larger than in the paramagnetic high spin state at 
room temperature^ This significant enhancement of the 
conductivity in the paramagnetic state agrees with Good- 
enough's hypothesis that thermal excitation populates 2 
orbitals which have a substantial degree of covalent mix- 
ing and itineracy, as opposed to the more localized A 
orbitals^ We see, however, that the magnetic suscep- 
tibility can be described by both a thermally induced 
Pauli susceptibility (Fig. 5) and a low spin to high spin 
transition^ Neutron scattering experiments would be 
very useful to quantify this and to correlate the phonon 
structure and frequency change with the thermally in- 
duced volume change. We conclude that the FeSb2 is 
a promising model system to study the applicability of 
the Kondo Insulator vs. the nearly itinerant magnetic 
semiconductor picture for 3d intermetallic compounds, 
in addition to FeSi. This is further supported by prelim- 
inary optical spectroscopy measurements^ 
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TABLE I: FeSb2 variable temperature refinements. Space group Pnnm, Fe atom at 2a site at (000), Sb atom at 49 she at 
(xyO), sample contains ~ 8wt.% of Sb impurity. Other values are available upon request 



T(K) 


20 


100 


150 


200 


250 


300 


400 


500 


a(A) 


5.82117(6) 


5.82333(11) 


5.82480(8) 


5.82626(8) 


5.82755(8) 


5.82764(9) 


5.82957(8) 


5.83298(7) 


b(l) 


6.50987(7) 


6.51336(12) 


6.51738(9) 


6.52304(9 


6.52911(9) 


6.53345(10) 


6.54221(9) 


6.55077(8) 


c(A) 


3.16707(4) 


3.17035(6) 


3.17525(5) 


3.18218(5) 


3.18977(5) 


3.19630(5) 


3.20785(5) 


3.21926(4) 


V(i 3 ) 


120.016(2) 


120.249(5) 


120.540(3) 


120.939(3) 


121.367(3) 


121.698(4) 


122.342(3) 


123.010(3) 


Sb(x) 


0.1875(1) 


0.1872(2) 


0.1875(2) 


0.1882(2) 


0.1882(2) 


0.1887(2) 


0.1894(2) 


0.1901(2) 


Sb(y) 


0.3554(1) 


0.3549(2) 


0.3550(2) 


0.3551(2) 


0.3558(2) 


0.3557(2) 


0.3564(2) 


0.3574(1) 


x 2 


6.740 


9.606 


8.822 


7.671 


8.842 


9.145 


8.186 


7.671 


Fe-Sb x 2 (A) 


2.5583(7) 


2.5558(9) 


2.5586(10) 


2.5628(10) 


2.5691(9) 


2.5708(10) 


2.5800(11) 


2.5904(9) 


Fe-Sb x 4 (A) 


2.5889(6) 


2.5931(8) 


2.5935(8) 


2.5934(8) 


2.5944(8) 


2.5951(9) 


2.5948(8) 


2.5945(7) 


Fe-Sb-Fe (o) 


128.95(2) 


128.99(3) 


128.99(3) 


128.98(3) 


128.83(3) 


128.86(3) 


128.74(3) 


128.58(3) 



TABLE II: FeSb2 high pressure refinements. Space group P„„ m , Fe atom at 2a site at (000), Sb atom at 4g site at (xyO), 
sample contains ~ 8wt.% of Sb impurity 



P (GPa) 


ambient 


0.79 


1.98 


2.91 


3.75 


4.15 


5.93 


7.14 


a (A) 


5.82788(9) 


5.8207(2) 


5.8042(2) 


5.7933(2) 


5.7843(2) 


5.7740(2) 


5.7583(2) 


5.7432(3) 


b(A) 


6.53336(9) 


6.5198(2) 


6.4907(2) 


6.4733(2) 


6.4599(2) 


6.4456(2) 


6.4217(2) 


6.3986(3) 


C (A) 


3.19626(5) 


3.1799(1) 


3.1570(1) 


3.1423(1) 


3.1312(2) 


3.1201(1) 


3.1012(1) 


3.0857(2) 


V (A3) 


121.699(4) 


120.676(6) 


118.933(6) 


117.843(5) 


117.001(6) 


116.118(4) 


114.677(5) 


113.394(6) 


Sb(x) 


0.1887(1) 


0.1858(4) 


0.1854(4) 


0.1857(3 


0.1839(4) 


0.1850(3) 


0.1842(4) 


0.1842(4) 


Sb(y) 


0.3561(1) 


0.3552(4) 


0.3537(3) 


0.3518(3) 


0.3527(3) 


0.3558(3) 


0.3585(4) 


0.3502(4) 


2 

X 


6.681 


2.625 


2.394 


6.126 


4.001 


19.91 


22.06 


12.89 


Fe-Sb x 2 (A) 


2.5733(7) 


2.556(2 


2.535(2) 


2.519(2) 


2.515(2) 
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